Superconducting critical field far above the Pauli limit in one-dimensional Lio 9M06O17 
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The upper critical field /fc2 of purple bronze Lio.gMogOn is found to exhibit a large anisotropy, in quantitative 
agreement with that expected from the observed electrical resistivity anisotropy. With the field aligned along 
the most conducting axis, Hc2 increases monotonically with decreasing temperature to a value five times larger 
than the estimated paramagnetic pair-breaking field. Theories for the enhancement of //^a invoking spin-orbit 
scattering or strong-coupling superconductivity are shown to be inadequate in explaining the observed behavior, 
suggesting that the pairing state in LIq gMogOn is unconventional and possibly spin-triplet. 
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Superconductivity in quasi-one-dimensional (qlD) conduc- 
tors has attracted sustained interest from the theoretical com- 
munity |1], largely due to the fact that under certain condi- 
tions, rare phenomena such as spin-triplet pairing |J2|-|4| or the 
spin-singlet, spatially inhomogeneous Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state [5-7] may be realized. The or- 
ganic conductors (TMTSF)2X (X = PFe, CIO4) have been 
most extensively studied in this regard, though the nature 
of their pairing state has not yet been fully determined. In 
(TMTSF)2PF6, the constant Knight shift across the supercon- 
ducting (SC) transition Tc [8], together with the observed vio- 
lation of the Pauli paramagnetic limit [9] supports triplet pair- 
ing, while in (TMTSF)2C104, a Kni ght shift suppression be- 
low Tc and the presence of line nodes [10] indicate collectively 
a £/-wave, spin-singlet pairing state. The possible realization 
of the FFLO state in (TMTSF)2C104 at low T, as suggested by 
recent angular studies of Hc2 1 1 1], is also consistent with sin- 
glet rather than triplet pairing. Theoretically, the coexistence 
of spin-(SDW) and charge-(CDW) density-wave instabilities 
can lead to a complex phase diagram where both singlet and 
triplet phases lie in close proximity, with triplet pairing be- 
coming dominant as repulsive interchain interactions are en- 
hanced isl I12I \1M ■ Intriguingly, the triplet state known to ex- 
ist in the q2D perovskite superconductor Sr2Ru04 might also 
arise from repulsive interactions between qlD bands lll4ll . 

Lio.gMogOn (LiMO) is a transition metal oxide with qlD 
electronic properties. It is metallic at high T, semiconducting 
below a temperature 15 K < Tmin < 30 K and superconduct- 
ing below Tc ^ 2K 11511 . Despite having a Tc higher than the 
(TMTSF)2X family, its SC properties have received little at- 
tention to date. While the presence of a density wave (DW) 
transition in LiMO was initially discounted, recent magneto- 
transport data appear to suggest some form of DW gapping 
III6II . The precise nature of the DW however, and its rela- 
tion to the superconductivity, has yet to be resolved. Finally, 
signatures of superconductivity have been found to emerge 
at high magnetic fields 11611 in LiMoO crystals that are non- 
superconducting in zero-field, consistent with theoretical pre- 
dictions for a qlD superconductor with triplet pairing |J6| tzl]. 

Here, we report a detailed temperature and orientational 
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FIG. 1: Zero-field resistivity curves for LiggMofjOn for l\\a,b and 
c, scaled by 2000 (Pc) and 100 (Pa) for clarity. Note that the PtiT) 
trace only goes down to 4.2 K, as explained in the Supplementary In- 
formation. Lower inset: Superconducting transition as seen in Pc{T) 
(blue) and its temperature derivative dpc/dT (brown). Upper inset: 
f-axis magnetoresistance Apc/pc on the same single crystal at T = 
4 K, i.e. just above Tc, plotted versus H^ (H||a). The slope gives a 
measure of the in-chain mean-free-path (see text). 



study of Hc2 in crystals that are superconducting in zero-field 
and extract Hc2{T) for fields applied along the three crystal- 
lographic axes. With the field parallel to the zigzag chains 
(Hll^j), Hc2 increases monotonically with decreasing temper- 
ature to a value five times larger than the usual Pauli param- 
agnetic limit. We show evidence that LiMoO is a strongly- 
coupled superconductor in the clean limit. However, the large 
Hc2 values can neither be explained wholly by the effects of 
strong coupling fl2'| nor by spin-orbit scattering as parame- 
terized by the Werthamer-Helfand-Hohenberg (WHH) theory 
11811 . Such a finding points to the possibility that the low-field 
SC state in LiMO might also be a spin-triplet. 
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FIG. 2: Field sweeps of the c-axis resistivity of Lio.gMogOn in ~ 0. 1 K steps forH aligned along the three crystallographic axes. 



Figure [T] presents zero-field resistivity p{T) curves from 
300 K down to 1.6 K (4.2 K) for I|la,c (l\\b) respectively, 
scaled to incorporate all three curves on the same set of axes. 
(Full details of the samples and measurement techniques are 
given in the Supplementary Information.) On cooling from 
room temperature, p(r) is metallic down to about 15 K, fol- 
lowed by a well-defined upturn and finally, superconductivity. 
The size of the resistivity upturn is much smaller than found 
in the (non-superconducting) crystals described in Ref. lll6ll . 
consistent with the anti-correlation between Tc and the size of 
the resistivity upturn first reported by Matsuda et al. \IV\. The 
lower inset shows a blow-up of the c-axis resistive transition at 
Tc, the mid-point of which, as defined using the maximum in 
dpf/dr, occurs at 2.2 K. The resistivity anisotropy extracted 
from this set of curves is 80:1:1600 (150:1:1600) for Pa'-9b'-9c 
at r = 300 K (4.2 K) respectively. This is almost two orders of 
magnitude larger than the anisotropics reported in the recent 
literature [20, 21 ], highlighting the extreme care needed to iso- 
late the individual current directions in such low-dimensional 
systems (see the SI section for further details). 

Figure |2] shows c-axis resistivity curves obtained in an 18 
Tesla pumped ^^He cryostat with the field applied along a, b, 
and c, for temperatures between 0.33 and 2.40 K. The Pc{H) 
curves for H\\b show an unusual broadening at intermediate 
temperatures, the origin of which is not understood at present. 
The large anisotropy in Hc2 is apparent from inspection of the 
field scales in the three different panels. Field alignment for 
H| I fe required an accuracy < 1° that was difficult to achieve in 
our ^He system and higher Hc2 values (H||/7) were observed 
in a second set of measurements performed on the same crys- 
tal in a pumped "^He system. The phase diagram for Hc2{T) 
obtained from fixed-field temperature sweeps in the latter is 
shown in Fig.|3] where the open squares, two-tone squares and 
open circles refer to measurements performed with H\\a,b,c 
respectively. Here, Hc2{T) is determined by the maximum in 
the relevant derivative dp^/dr or dpf/d//. As shown in the SI 
section, choosing a different criterion does not change qualita- 



tively the overall behavior, nor the anisotropy parameters . All 
data in Figs. |2] and [3] were corrected for the remnant field by 
symmetrizing with respect to positive and negative field val- 
ues. Data for a second single crystal with a slightly lower Tc 
and Hc2{0) are shown (for H\\b only) in the SI section. 

From the initial slopes of the increase in Hc2{T) below Tc 
(dashed lines in Fig. O, we obtain dHc2l<^T\T<T^ values of - 
1.5 T/K, -19.5 T/K and -0.5 T/K and a corresponding critical 
field ratio of 3:39:1 for H||fl:fe:c respectively. According to 
anisotropic Ginzburg-Landau (GL) theory 
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where subscripts ij refer to crystalline axes, superscripts 
i,i refer to the direction of the applied magnetic field, ^,j 
(o= \'i,j) is the orientation dependent coherence length, Oa/ajj 
(o= vjj/vj:) is the anisotropy in the diagonal elements of the 
conductivity tensor and vtj are the respective Fermi veloci- 
ties. Squared, the SC anisotropy is approximately 9:1500:1, 
or 170:1:1500 when inverted. Note that these ratios are in 
excellent quantitative agreement with those (150:1:1600) ob- 
tained from the normal state resistivity measurements. 

A more complete Hc2 (T) phase diagram, extracted from the 
data plotted in Fig. |2] is shown in Fig. |4] (Note that the Hc2 
values for H||fl,c have been re-scaled in this plot). Using the 
(linearly) extrapolated zero-temperature values Hc2{0) from 
the phase diagram and the equation 
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we obtain estimates for the three coherence lengths; ^^(0) ~ 
300 A, t,JO) ~ 100 A and ^,.(0) ~ 25 A. Significantly, the 
interchain coherence lengths are both longer than the lattice 
spacing (or more precisely 2^a(^.)(0) > a{c) — d 11221] . where 
a = I2.73A and c = 9.51 A are the a, c-axis lattice constants 
and d ^ 3 A is the approximate width of the M0O4 octahedra). 
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FIG. 3: Phase diagram for Hc2 as a function of temperature for fields 
applied along the three crystallographic directions and temperatures 
down to 1.6 K obtained with the He vapor pressure system. 



implying the absence of Josephson coupling and a continuous 
phase of the SC order parameter across the chains. Hence, de- 
spite the extreme one-dimensionality of LiMO in the normal 
state, its superconductivity appears to be described satisfacto- 
rily using anisotropic-3D GL theory. 

For a qlD spin-singlet superconductor with H|| chain, or- 
bital pair-breaking is minimized (due to the small interchain 
electron velocities) and superconductivity can only be de- 
stroyed once the Zeeman energy arising from spin-splitting 
of Cooper pairs exceeds the SC condensation energy. For an 
isotropic BCS superconductor, the Pauli limit is expressed as 
fJoHp = 1.847;: ~ 4.0 T for Tc = 2.2 K. Hp can also be cal- 
culated independently 12311 using actual values (see SI sec- 
tion) for the Pauli susceptibility Xp (= 2.8 x 10^^) and for 
the condensation energy Uc (= 2.2 mJ/mol, estimated from 
the specific heat anomaly at 7^) obtained on crystals taken 
from the same batch and with similar 7^. values. These give 
/JoHp ~ 3. 1 T for LiMO, i.e. comparable with the BCS value 
but still five times smaller than the measured 7/c2(0) for ll\\b. 



According to WHH theory OlSll . spin-orbit scattering (e.g. 
at the Mo site) can act to limit Zeeman splitting and thus in- 
crease the value of Hc2 beyond the usual Pauli limit. It is 
expressed using two dimensionless parameters, the Maki pa- 
rameter a and the spin-orbit scattering A-so- The former is 
constrained through the expression a ~ 0.528 dHi.2/dT\j'<j'^ 
(= 7.5 for the ^He data), while A-so = 2h/3nkBTcTso (= 32Ms 
determined from the best curve fit to the data, shown in Fig. |4] 
by a dashed line. The associated spin orbit scattering time Tso 
can be converted to a mean-free-path ^so using the measured 
(in-chain) Fermi velocity [24']. For both data sets (i.e. from 
the ^He (Fig.H and ^He cryostats (Fig.©), £so - 120 A. 

An estimate for the transport mean-free-path £q can be ob- 
tained from the low-T interchain magnetoresistance (MR). 
With I||c and H||fl, Boltzmann theory gives for a qlD metal in 
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FIG. 4: Phase diagram for //c2(7') for temperatures down to 0.3 K 
obtained with the He cryostat. The data were obtained from the field 
derivative of the magnetic field sweep curves shown in Fig.|2] Note 
that for H||a,c, the H1.2 values have been normalized to emphasize 
the similarities of their T-dependencies. The dashed line is a fit to 
the WHH theory with a = 7.5 and A-so - 32 (see text for details). 
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As shown in the top inset of Fig. [T] Apc/pc °= B^ at low 
fields with a slope that yields Iq ~ 650 A at T = 4 K. Ac- 
cording to the Abiikosov formula (1251, Xq/Tso = (Z/137)^ (Z 
= 42 for Molydenum). Thus, one expects Tso to be ~ 100 
times longer than the transport lifetime Tq, in marked con- 
trast to what is obtained from the WHH parameters. More- 
over, the diffusion constant obtained from the Maki parameter 
gives io < 1 A, showing clearly that WHH theory is inapplica- 
ble here (it is only truly valid in the dirty limit). It also shows 
that the spin-orbit scattering alone cannot account for the val- 
ues of Hc2 observed in LiMO. 

Strong (e.g. electron-phonon) coupling can also act to re- 
duce the effects of Pauli limiting through renormalization of 
the band splitting. Fits to the specific-heat anomaly yield a 
coupling constant A = 1.2 ± 0.1, intermediate between those 
of Nb and Pb. Even with this strength of coupling however, the 
magnitude of //c2(0) in LiMO is more than double the renor- 
malized value (see SI section for more details) II12II . While 
one cannot exclude the possibility that the high Hc2 found 
in LiMO is due to a combination of singlet pairing, strong- 
coupling superconductivity and strong spin-orbit scattering, it 
would require a precise combination of aU of these effects to 
realize the present situation. 

In a triplet superconductor, the spins of the Cooper pair can 
be co-aligned, making the Pauli pair-breaking effect redun- 
dant. Triplet superconductivity is, however, extremely frag- 
ile and can easily be destroyed by impuiities ll26ll . For our 
samples, £0 > 2^^,(0), i.e. within the clean limit, and while 
the FFLO state cannot be ruled out at the lowest temperatures 



and highest fields, significantly we find in LiMO, three of the 
key ingredients for realizing qlD triplet superconductivity: 
extreme one-dimensionality, a H^iiT) profile with minimal 
paramagnetic limiting and a long mean-free-path. The recent 
observation of scaling in the longitudinal fe-axis MR in LiMO 
111611 provides compelling evidence that some form of DW de- 
velops below Tmin. As discussed in the introduction, the nature 
of the DW fluctuations near Tc may ultimately determine the 
pairing state. In qlD systems close to a Peierls-type CDW 
instability, i-wave is the dominant pairing channel [27]. How- 
ever structural [.28,1 . thermodynamic [,1 1.[ and optical studies 
II6D have all failed to find evidence of a genuine phase tran- 
sition in LiMO at T = T^m- We also find no evidence of 
a specific heat anomaly at T^m, contrary to an earlier report 
03011 . This lack of evidence has lent support to the notion that 
the CDW instability in LiMO is in fact driven by electronic- 
interactions [28[, as is the case for strongly-interacting cou- 
pled Luttinger liquids 113 1|[ . with either singlet or triplet pair- 
ing competing for the ground state. 

Finally, for H±fe, Hc2{T) displays pronounced upward cur- 
vature, leading to a reduction in the SC anisotropy with de- 
creasing temperature. A similar enhancement has been ob- 
served in (TMTSF)2PF6 and attributed to the formation of in- 
sulating (SDW) and SC domains at pressures approaching the 
SDW phase 13211 . This model, in which Hc2 is determined 
by the largest penetration depth perpendicular to the applied 
field, can account for the very similar Hc2 values found for 
H|lfl and H||fe in qlD (TMTSF)2PF6 as well as the upturn in 
Hc2 for H||c. In LiMO however, there is no clear evidence for 
such an SDW phase at these temperatures and field scales and 
the anisotropy in Hc2, though reduced with decreasing tem- 
perature, always remains large. 

An alternative explanation for the enhancement of Hc2 is 
a field-induced reduction in the effective dimensionality of 
the electronic system, as observed in other qlD conductors 
llllll33M35[[ . ForH_Lfe, the semiclassical motion of quasiparti- 
cles along, say, k oscillates with an amplitude tii/evpB, where 
tk is the ^-axis hopping parameter (Note that even though 
Pc{T) is non-metallic at low-T, the positive, quadratic trans- 
verse MR shown in the inset to Fig. [T] indicates that coher- 
ent quasiparticles do exist in this temperature regime.) As 
B increases, the amplitude of the oscillatory interchain mo- 
tion decreases [36]. This gradual confinement suppresses or- 
bital pair-breaking and leads to an overall decrease in v^^^) and 
in turn t,^, which now becomes field-dependent. As a result, 
H'^.2 increases with decreasing temperature. While this field- 
induced dimensional crossover will occur in both orthogonal 
field orientations, it is expected to happen at different field 
scales [37]. Moreover, it remains an open question whether 
this confinement process can lead to an enhancement of Hc2 
at fields well below the crossover field. In the case of triplet 
pairing, this field-induced dimensional crossover can lead ulti- 
mately to re-entrant superconductivity [6, 7], though in LiMO, 
such behavior has only been seen to date with the field aligned 
parallel to the molybdate chains and only in samples that do 
not superconduct in zero-field lll6ll . In the present batch of 



crystals, we observe a strong negative magnetoresistance for 
tl\\b, but as yet, no sign of re-entrant superconductivity. 

In conclusion, we have demonstrated the existence of a 
highly anisotropic, yet still three-dimensional SC state in 
LiMO with an anisotropy (near Tc) that is in excellent quan- 
titative agreement with the measured normal-state electrical 
anisotropy. The magnitude of the b-axis upper critical field 
exceeds the usual Pauli limit by a factor of five. We have 
shown that neither spin-orbit scattering nor strong-coupling 
superconductivity can account fully for this enhanced Hc2 and 
thus have speculated that LiMO is a viable candidate for the 
realization of triplet superconductivity. Although the effective 
dimensionality of the electronic state just above Tc is yet to 
be determined, LiMO displays all the hallmarks of a ID Lut- 
tinger liquid at least above T = Tj^in 12413811 . suggesting that 
the superconductivity may in fact have a higher dimensional- 
ity than the normal state out of which it condenses. In this 
case, pairing has to involve electrons on different chains, thus 
providing a test-bed for theoretical claims that triplet pairing 
in qlD superconductors is stabilized in the presence of (repul- 
sive) interchain interactions [|3i ll2lll3ll . 

We thank J. F Annett, P Chudzinski, B. Gyorffy, R. H. 
McKenzie, J. Merino and N. Shannon for stimulating discus- 
sions. This work was supported by the EPSRC (UK). XX 
acknowledges financial support from NSFC (No. 1 1 104051), 
NEH a Royal Society Wolfson Research Merit Award. 



[1] V. J. Emei-y, Synth. Met. 13 21 (1986). 

[2] D. Podolsky, E. Altman, T. Rostunov and E. Demler, Phys. Rev. 

Lett. 93, 246402 (2004). 
[3] J. C. Nickel, R. Duprat, C. Bourbonnais and N. Dupuis, Phys. 

Rev. Lett. 95, 247001 (2005). 
[4] K. Kuroki and Y. Tanaka, J. Phys. Soc. Japan 74, 1694 (2005). 
[5] A. I. Buzdin and V. V. Tugushev, Zh. Eksp. Teor. Fiz. 85, 735 

(1983) [Sov. Phys. JETP 58, 428 (1983)]. 
[6] A. G. Lebed, Pisma Zh. Eksp. Teor. Fiz. 44, 89 (1986) [JETP 

Lett. 44, 114(1986)]. 
[7] N. Dupuis, G. Montambaux, and C. A. R. Sa de Melo, Phys. 

Rev. Lett. 70, 2613 (1993). 
[8] I. J. Lee etai, Phys. Rev. Lett. 88, 017004 (2002). 
[9] I. J. Lee, P M. Chaikin and M. J. Naughton, Phys. Rev. B 62, 
R14669 (2000). 
[10] J. Shinagawa et al., Phys. Rev. Lett. 98, 147002 (2007). 
[11] S. Yonezawa et al., Phys. Rev. Lett. 100, 1 17002 (2008). 
[12] K. Kuroki, R. Arita and H. Aoki, Phys. Rev. B 63, 094509 

(2001). 
[13] Y. Fuseya and Y Suzumura, J. Phys. Soc. Japan 74, 1263 

(2005). 
[14] S. Raghu, A. Kapitulnik and S. A. Kivelson, Phys. Rev. Lett. 

105, 136401 (2010). 
[15] M. Greenblatt et al.. Solid State Commun. 51, 671 (1984). 
[16] Xiaofeng Xu etai, Phys. Rev. Lett. 102, 206602 (2009). 
[17] J. P Carbotte, Rev. Mod. Phys. 62, 1627 (1990). 
[18] N. R. Werthamer, E. Helfand and P C. Hohenberg, Phys. Rev. 

147, 295 (1966). 
[19] Y Matsuda, M. Sato, M. Onoda and K. Nakao, J. Phys. C 19, 
6039 (1986). 



[20] M. S. da Luz et al, Phys. Rev. B 76, 233105 (2007). 

[21] H. Chen et al., Europhys. Lett. 89 67010 (2010). 

[22] N. E. Hussey et ai, Phys. Rev. B 53, 6752 (1996). 

[23] E Zuo et ai, Phys. Rev B 61, 750 (2000). 

[24] E Wang et ai, Phys. Rev Lett. 96, 196403 (2006). 

[25] A. A. Abrikosov and L. P. Gor'kov, Zh. Eksp. Teor. Eiz. 42, 

1088 (1962) [Sov Phys. JETP 15, 752 (1962)]. 
[26] A. P Mackenzie et ai, Phys. Rev Lett. 80, 161 (1998). 
[27] H. Bakrim and C. Bourbonnais, Europhys. Lett. 90, 27001 

(2010). 
[28] C. A. M. dos Santos et al., Phys. Rev Lett. 98, 266405 (2007). 
[29] J. Choi et al., Phys. Rev B 69, 085120 (2004). 
[30] C. Schlenker, H. Schwenk, C. Escribe-Eilippini and J. Marcus, 

Physical35B, 511 (1985). 
[31] D. Boies, C. Bourbonnais and A.-M. S. Tremblay, Phys. Rev. 

Lett. 74, 968 (1995). 
[32] L J. Lee, P. M. Chaikin and M. J. Naughton, Phys. Rev Lett. 

88, 207002 (2002). 
[33] K. Behniaef a/., Phys. Rev Lett. 74, 5272 (1995). 
[34] L J. Lee, M. J. Naughton, G. M. Danner and P. M. Chaikin, 

Phys. Rev. Lett. 78, 3555 (1997). 
[35] N. E. Hussey et ai, Phys. Rev Lett. 89, 086601 (2002). 
[36] L. P Gor'kov and A. G. Lebed, J. Phys. (Paris) Lett. 45, 433 

(1984). 
[37] A. Narduzzo et al.. New J. Phys. 8, 172 (2006). 
[38] N. Wakeham et ai, Nat. Commun. 2, 396 (201 1). 



SAMPLE CHARACTERIZATION AND RESISTIVITY 
MEASUREMENTS 

Single crystals of Lio.glVloeOiy were grown using a tem- 
perature gradient flux method [1]. The crystals reported in 
this Letter were all selected from the same batch and face- 
indexed using a single-crystal X-ray diffractometer to check 
their crystallinity and to identify the a- and b-axes. The lattice 
parameters were a = 12.73 (± 0.01) A, b = 5.53 (± 0.01) Aand 
c = 9.51 (± 0.01) A, while the monoclinic angle (p) was de- 
termined to be 90.6°, in agreement with previous crystallo- 
graphic studies 021 13|]. Once face-indexed, the selected sam- 
ples were then cut into rectangles of appropriate dimensions 
along the a- and b-axes and cleaved to obtain clean surfaces 
onto which electrical contacts were made. 

For the Pa(Pc) measurements, we used the so-called quasi- 
Montgomery method in which the resistivity is measured in 
only one (the most resistive) configuration. In both cases, the 
dimensions of the crystal were chosen such that the a- or c- 
axis dimension of the equivalent isotropic crystal, estimated 
by multiplying the true a- or c-axis dimension by the square 
of the resistive anisotropy [4], was always at least one order of 
magnitude longer than the corresponding length along b. This 
ensured that the current flow was uniaxial along a or c. 

For the Pa{T) measurements, crystals were cut into squares 
of approximate dimensions 300 x 300 x 20 /jm^ and four gold 
wires were attached with Ag paint at the corners of each crys- 
tal, which itself was raised above the sample platform with 
the current/voltage contacts acting as posts in order to short 
out the highly resistive c-axis component. For the Pc{T) mea- 
surements, thicker samples were employed, typically 200 x 
200 X 50 ijnv' with the gold wires painted onto the faces of the 
crystal in such a way as to short out the more resistive a-axis. 
In both cases, the interchain resistivities were measured using 
ac lock-in detection techniques with excitation currents of 1 
- 10 /jA. No non-linear effects were observed, neither in the 
normal state nor in the mixed state, over this range of currents. 

Of all the individual components of the electrical conduc- 
tivity tensor, the most difficult one to measure reliably is the 
in-chain conductivity (resistivity) since any spurious voltage 
drop perpendicular to the chain direction can lead to widely er- 
roneous resistivity values and an incorrect temperature depen- 
dence. While the Montgomery technique is a reliable method 
for determining the inter-chain components, in our experi- 
ence, it is not an effective means of isolating the small in-chain 
resistivity since its applicability relies on the ability to short 
out the voltage drop in the third orthogonal direction which 
often cannot be guaranteed. In our experiments therefore, we 
choose to mount our electrical contacts in such a way that 
the voltage drop in both orthogonal directions is shorted out 
rendering the current flow between voltage contacts strictly 
one-dimensional. This is achieved by having current pads that 
completely cover the ends of the (bar-shaped) sample, using 
either conductive paint or evaporated gold strips, and voltage 
pads that are painted (or evaporated) across the entire width 
of the sample as well as down the sides. A detailed descrip- 
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FIG. 5: Demonstration of the different criteria used to determine Hc2 
in Lio.gMosOn. Tlie solid blue line is Pc{H) field-sweep data ob- 
tained at r = 0.81 K with H\\b. The 10% and 50% criteria refer 
to the field scale at which the resistivity reaches 10% and 50% re- 
spectively of the maximum in the Pc{tl) curve. The crossing-point 
criterion corresponds to the field scale at which extrapolations of the 
mixed-state and normal-state resistivities meet. Note that the 50% 
criterion coincides with the maximum in dpc/dT, as expected. 



tion of our determination of Ph{T) has already been published 
elsewhere, and we direct the interested reader to the Supple- 
mentary Information section of Ref. [5] for further details. 
Typical dimensions of our Pb{T) crystals were 750 x 150 x 20 
jum^ and the excitation currents in this case ranged between 
lOandlOO^^A. 

We note that previous resistivity data reported in the liter- 
ature and performed using the Montgomery method have a 
much lower electrical anisotropy than ours, with ph{T) dis- 
playing a sub-linear temperature dependence at high T remi- 
niscent of the a- or c-axis resistivities. Of course, a low resis- 
tivity value (for current along the chains) does not automati- 
cally guarantee that one has isolated the in-chain current re- 
sponse, since current paths can sometimes be such that lower 
(and even negative) resistivity values are obtained artificially. 
However, a strong indicator that we are measuring the intrinsic 
b-axis response is the reproducibility of our data for a number 
of different samples. Again, we refer the interested reader to 
the Supplementary Information of Ref. |,5J for more detail. In 
total, we obtained Pb{T) curves for over 30 single crystals to 
allow us to better identify the intrinsic T-dependence of the 
in-chain resistivity and to determine its absolute magnitude. 
It is important to acknowledge that the resistivity anisotropy 
we obtain is not only in good agreement with measured Hc2 
anisotropy, but also with that measured independently by op- 



tical conductivity H. 

Finally, note that the pb{T) curve displayed in Figure 1 in 
the main manuscript and the one we believe to be the most 
representative of the intrinsic b-axis resistivity, has a lower 
limit of only 4.2 K. This is because the data were taken in a 
^He system that could not be pumped. Within a single batch 
however, the 71- values (and the size of their resistive upturns) 
tend to vary by less than 10% (40%) respectively. Moreover, 
as shown in a later section of the Supplementary Information, 
we observed a large specific heat anomaly at Tc = 2.2 K in 
one of these crystals, confirming that the superconductivity is 
indeed bulk and three-dimensional. 



DETERMINATION OF Hc2 



In Figures 3 and 4 of the Letter, we determine Hc2 using 
the maximum in the temperature or field derivative of Pc- In 
this section, we consider the effect of adopting different cri- 
teria in the determination of the superconducting anisotropy 
in Lio.gMogOn. Figure SI shows a typical pc{H) sweep for 
H||fe. The temperature of the sweep is T = 0.81 K. As might 
be expected, the maximum in dpc/dT gives a near-identical 
Hc2 value to that determined by the midpoint of the transition. 
The 10% criterion, i.e. the field at which p^ reaches 10% of its 
maximum value, is relatively straightforward to identify and 
can be used to set a lower limit on Hc2{T). For the upper 
limit, the 90% criterion is much harder to define due to the 
broad fluctuation regime (for M\\b). However, if instead we 
use the crossing point between the mixed-state and normal- 
state resistivities for the definition of our high Hc2 value (the 
so-called 'Junction' criterion of Lee et al. |7|]), as illustrated 
in Figure SI, we obtain the three sets of curves shown in Fig- 
ure S2. Consistent data were obtained on a second crystal, as 
described in the following section. 

As can be seen, the T-dependence ofHc2 is only slightly de- 
pendent on the criterion used. More importantly, we find that 
the anisotropy of the initial slopes is relatively insensitive to 
the different criteria. After scaling the Hc2 values obtained for 
H||fo in the ^He system to those obtained in the pumped ^He 
system (due to a slight misalignment of the applied field, as 
explained in the main text), we obtain H||a:fe:c dHc2/dT\T<Tc 
values of -l.l:-17:-0.4 T/K and -2.4:-30:-0.7 T/K for the 10% 
and crossing-point criteria respectively, corresponding to crit- 
ical field ratios of 2.8:44:1 and 3.5:44:1 for H||a:fe:c respec- 
tively. These compare favorably with the ratio 3:39: 1 obtained 
using the maximum in dpc/dT, affirming that the anisotropy 
ratios are robust. Finally the linear extrapolations of the differ- 
ent Hc2{T) curves plotted in Figure S2 provide the following 
estimates of the uncertainty in our Hc2(0) values: H"2{0) = 5.0 
± 0.5 T, ///:;, (0) = 15.2 ± 0.2 T, H;:.2{0) = 1.2 ± 0.2 T. 
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FIG. 7: The temperature dependence of the upper critical field 
HciiJ) for a second crystal {Tc = 2.1 K) and the same crystal de- 
scribed in the main manuscript. 



Hc2 OF A SINGLE CRYSTAL WITH A LOWER T,- 

In our continuing search for field-induced superconductiv- 
ity in this new batch of crystals, we found similar Hc2{0) val- 
ues (i.e. > 13 Tesla) in at least two other samples. Figure S3 



shows the results of a more detailed study of the Hc2{T) pro- 
file for H||fe for a second sample with a slightly lower Tc (= 
2.1 K) and a correspondingly lower Hc2{0) (=13 T). Super- 
imposed on this curve are data for the same crystal described 
in the main manuscript (these are the same data as shown in 
Figure 4), scaled by a factor of 0.92. Thus, for a 5% reduction 
in Tc, there is a corresponding 8% reduction in Hc2{0). The T- 
dependence of the upper critical field however is essentially 
identical in the two samples. 



ESTIMATING THE PAULI LIMIT FROM HEAT CAPACITY 
AND MAGNETIC SUSCEPTIBILITY MEASUREMENTS 

Figure S4 shows heat capacity data, divided by tempera- 
ture, for a single crystal from the same batch as the crystal 
described in the main text. Measurements were conducted us- 
ing a long-relaxation calorimeter [10] on a sample with mass 
~ 4.4 mg. A large specific heat anomaly, associated with 
the superconducting transition, is observed centered at T = 
2.15 K, in excellent agreement with the resistive transition re- 
ported in Figure 1. The size of the jump in C at Tc, is far in 
excess of that expected for weak-coupling superconductivity 
(AC/'YTc=2.07), so we have fitted the data using the strong- 
coupling alpha model [18]. Here the temperature dependence 
of the energy gap is assumed to be the same as in the weak- 
coupling s-wave BCS model, but the absolute value is mul- 
tiplied by a constant a to account for the enhanced value of 
A/Tc due to strong coupling. This model has been found to 
account well for data, both for conventional strong-coupling 
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FIG. 8: The low-temperature specific heat of a Lio.gMogOn single 
crystal, y, divided by temperature. The solid line is a fit to the alpha- 
model for a strongly-coupled superconductor as described in the text. 
The dashed line is a guide to the eye. 



superconductors fS], as well as more exotic materials such as 
MgB2 [9] and the organic superconductors |10|. The fit gives 
the only free parameter in this model, a = 1.6. By integrat- 
ing the fitted curve, we get a superconducting condensation 
energy of Uc = 2.2 mJ/mol (^11 J/m^). In many different 
electron-phonon coupled superconductors there is a universal 
relation between AC/yTc and the ratio of Tc to the average 
phonon-frequency (Oin, which is approximately given by lll2ll 



AC 



1.53 



1+53 



coin 



In 



Win 

3r, 



This gives, COin/Tf. = 0.075, which is intermittent between Nb 
and Pb. Using this we can solve the McMillan equation 



1.271= ©inexp 



-1.04(1+ A.) 
X-IJ*{l+0.62X) 



to obtain a value for the coupling constant A, = 1.23 for /j* ~ 
0.15(A,= 1.12for^*=0.12). 

The inset to Figure S5 shows raw magnetic susceptibility 
X{T) data for another single crystal from the same batch mea- 
sured between 2 K and 200 K in a magnetic field of 5 Tesla 
applied along the b-axis. The same crystal also shows a su- 
perconducting transition in lower fields at Tc = 2.0 K. Plot- 
ted alongside this data set are earlier x(^) d^t^ reported by 
Matsuda et al. back in 1986 [11]. While the high temper- 
ature values are similar, the data of Matsuda et al. show a 
much larger enhancement of low T. This enhancement was 
attributed to an extrinsic Curie paramagnetic term caused by 
a concentration of 0.15 x 10^^ localized (S = 1/2) moments 



FIG. 9: Main panel: Residual magnetic susceptibility data for 
Lio,9Mo50i7 obtained from subtraction of the l/T Curie term. 
The dashed line is a guide to the eye. Inset: Raw susceptibility 
data for our Lio.gMosOn single crystal (solid symbols) and for a 
Lio.gMosOn single crystal reported by Matsuda et al. 1 1 11 (open 
circles). 



per conduction electron. In our crystals, the enhancement is 
a factor of 20 smaller Subtracting this small l/T term from 
the raw data, we obtain the curve plotted in the main panel of 
Figure S5. The residual susceptibility data show a very weak 
T-dependence, extrapolating to a zero-temperature value x(0) 
= 2.8 X 10^^. The ionic diamagnetism Xciia and the orbital van 
Vleck contribution 5(„. to x{T) are believed to be comparable 
in LingMofiOn, with x^/m having the slightly larger magni- 
tude fll]. This implies that the residual susceptibility data is 
a lower estimate of the full Pauli paramagnetic susceptibility 
term X/>- We can also obtain an estimate for Xp from the mea- 
sured Sommerfeld constant y(0) (= 1.6 mJ/mol.K^), assuming 
a Wilson ratio Rw = 2 (that is consistent with most strongly 
correlated metals as well as Luttinger liquids with repulsive 
interactions). From these considerations, we obtain Xp = 3.0 
X 10^'', in good agreement with the measured value. 

Combining our estimates for Uc and %/>, we obtain an esti- 
mate for the Pauli limit of Bp = (l/JoUc/Xpf-^ = 3.1 Tesla, i.e. 
lower than the limit estimated simply from the value of Tc. 
This estimate can be revised upward however if one consid- 
ers the effects of strong coupling. According to Figure 82 of 
Carbotte's review of strong-coupling superconductivity lll2ll . 
given the slope of HciiJ) near Tc (~ 14 T/K, as measured in 
the ^He system) and our estimate of A,, one expects a corre- 
sponding Pauli limiting field of 7.3 Tesla in the clean limit. 
However, the actual value of HdiO) in our crystals (> 15 T) 
is more than a factor of 2 larger than this strong-coupling ex- 
pectation. 
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